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Abstract 

We discuss a very peculiar subgroup of gamma-ray bursts among the BATSE sources. These 
bursts are very short (Tgo <0.1 s), hard, and came predominantly from a restricted direction of the 
sky (close to the Galactic anti-center). We analyze their arrival times and possible correlations, as 
well as the profiles of individual bursts. We find no peculiarities in the arrival times of Very Short 
Bursts (VSBs) despite their highly non-uniform spatial distribution. There is no dependence in 
the burst shapes on location. Bursts coming both from the burst-enhancement Galactic Anticenter 
region and from all other directions show considerable dispersion in their rise and fall times. 
Significant fraction of VSBs have multiple peaks despite their extremely short duration. Burst 
time properties are most likely to be consistent with two origin mechanisms: either with binary 
NS-NS mergers with low total masses passing through a phase of hypermassive neutron star, or 
with evaporation of the primordial black holes in the scenario of no photosphere formation. 
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1. Introduction 

Observations of gamma ray bursts by BATSE satellite allowed to accumulate the data set 
large enough for very detailed statistical studies. The distribution of the burst duration, defined 
as Tgo (time-interval during which 90 p e r cent of the fluen c e is a ccumulated), shows clearly 



bi-modal distribution (IKouveliotou et all [1993; Norri s et al. , 2000). Long bursts (Tgo > 2 s) 



show usually rather complex time-profiles, their energy spectra are generally softer, and a num- 
ber of these bursts are now identified with the X-ray, optical and radio counterparts. Studies 
of the detected a fterglows ( s ee J. Greiner's web pagaM) clea r ly support the cosmological origin 



of these bursts (Costa et all 1 1997c Ivan P aradiis et all 1 1997c iMetzger et all Il997l) . The hyper 



nova scenario, advocated strongly by Paczvnski ( 1995), is the most plausible explanation of this 
phenomenon (for a review, see Piran , 2000|: [Meszaros], 20021). For shor t bursts, the most pop- 
ular scenario is compact obj ect mergers dLattimer and Schramm . 19761) : for a review see e.g. 
Bulik and Belczvnski I (120041) . 

Recently, Zha ng et al.l (|2009) discussed the dichotomy of the two GRB populations and in- 
troduced the nomenclature of Type I/II bursts, instead of the short/long bursts. This was to denote 



1 http://www.mpe.mpg.de/ jcg/grbgen.html 
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the physically distinct categories of cosmological GRBs and to stress that the duration and hard- 
ness of the GRBs do not have to be the reliable indicators of their true physical nature. 

However, the exis tence of a s epar ate third group of extremely short bursts of still another 



origin is also likely. Cli ne et al.l dl999) suggested that very short bursts, lasting below 0.1 s, 



form a separate special subgroup of short bursts. They argued that the very short bursts have 
euclidean spatial distribution, rather different from the spatial distribution of bursts lasting 0. 1 
- 2 s (as studied later e.g. bv iGuetta and Piran" . 20051) . They suggested a local origin of these 
shortest bursts, possibl y due to the evap oration of primordial black holes. The argument was 
further strengthened by Clin e et al.l (12003b . and by ICline et al. (2006) where they stress excep- 
tional hardness of those events at the basis of KONUS data. These bursts were also found to have 



clearly anisotropic distribution across the sky (ICline et all 12001) 



Careful analysis of the spatial distribution of short bursts (Tgo < 2 s) was performed by 
iMagliocchetti et aL ( 2003 ). They showed an enhanced correlation between the short bursts at 
angular scale of 2-4°, which might suggest that up to ~ 1 3 % of the bursts may represent repeated 
activity. The anisotropy in the burst d istribution in a form of a spot around Galactic coordinates 
1 = 1 15° and b=30° was argued for bv lLitvin et all d200ll) . 

The properties of the short burst class are reviewed e.g. bv iNakar I (120071) . The classical 
distinction between the short and long bursts is defined by the Tgo < 2 s, and their duration 
distribution peaks at Tgo as 0.8 s. The temporal structure of these bursts is difficult for studies, 
because of the low signal-to-noise ratio and limited time resolution. Therefore, any more detailed 
studies of the temporal structure of the short bursts, performed in the case of the long GRBs, were 
not c arried out so far , becau se of the difficulty in analysis of their lightcurves. 

Nakar and Piran I (120021) find that most of the bursts from their sample of 33 bright events 



exhibit a rapid time variability on timescales much shorter than their total durations. In their 
sample, more than a half of bursts show at least two well separated pulses, and more than a third 
show rapid variability. For the single pulses, their durations range from 5 to 300 ms, with a peak 
around 50 ms. 



The pulse properties of 55 short bursts from BATSE were also studied by McB reen et al. 



(2001) who found that the rise and fall times as well as the FWHM and amplitudes show lognor- 
mal distributions. The median value of the rise time of a single pulse given by these authors is 
0.035 s, while for the fall time it is 0.056 and the asymmetry ratio is fese/fan = 0.65. The median 
total duration time for this sample was 0.095 and the median number of pulses was 2.5. 

A possible sub-class of t he short gamma ray bursts are the soft gamma-ray repeaters (SGR: 
Woods and Thomson , 2006). Their giant flares of durations about 0.1 s and hard (~ 300 keV) 



energy radiation are followed by a long (about 300 s), pulsating tail due to the rotation of a 
magnetized neutron star. The latter is much less luminous, and might eas ily be missed by the 
BATSE detector, in which case the event would be classified as a short GRB. lTanvir et al. Id2005l) 
found a correlation between the locations of short GRBs and positions of galaxies in the local 
Universe, indicating that 10-25% of the events originate locally, at z < 0.025. However, the 
signal was st ronges t for T gp > 0.5 s, which is longer than a typical duration of a SGR spike. 
Furthermore. lOfek I J2007I) finds, that the upper limit on the fraction of SGRs out of the BATSE 
sample is /sgr < 0.14. 

In the present paper we revisit the problem of the very short bursts (VSB), with the aim to 
test their timelike characteristics. We use the BATSE sample since the number of VSB with 
known localization on the sky in SWIFT data (only 8 events) is not high enough for our study. In 
Section|2] we present our sample of very short bursts from BATSE. We study their properties such 
as spatial distribution, arrival times and possible correlations, as well as the profiles of individual 
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bursts. In Section [3] we discuss the possible origins of the special properties of our subsample 
and the most plausible mechanism of emission. In Section|4]we conclude our work. 

2. Properties of the VSB from BATSE 

Properties of the bursts detected by BATSE are provided by the public current catalog It 
contains 2702 events. Out of those, 2040 bursts have the measurements of Tw as well as the 
coordinates in the sky. 

The subsample of short bursts (hereafter SB), with Tgo < 2, consists of 499 bursts. Out of 
those bursts, 51 have 790 < 0.1 s, and they further form our subsample of VSB. As shown in 
Cline et al. (2003), 21 of those come from 1/8 of the sky, extending from +90 to + 180 deg in 
galactic longitude and from -30 to + 30 deg in the galactic latitude (roughly in the direction of 
the Galactic Anticenter) so their spacial distribution is highly non-uniform. We now study other 
properties of this subsample which may make them different from the usual bursts and indicate 
their nature. 

2.1. Arrival times 

We construct a histogram of the arrival time distribution of VSB at yearly basis. The result 
is shown in Fig. Q] The data for all years were used, so in order to include as many bursts 
as possible we count the years from the beginning of the BATSE operation, i.e. the first year 
extends from May 1991 till April 1992 etc. There seems to be a drop in the number of VSB after 
1997, and the distribution peaks around the years 1993 - 1995. However, the probability that the 
histogram describes the constant rate can be accepted at a 0.89 confidence level using the KS test 
although the straight line representing the mean rate (5.7/year) gives a poor fit to this distribution 
{ X 2 = 13.5for8d.o.f.). 

Any departure from a constant rate would be rather surprising effect. Since the bursts seem 
to be grouped, we also construct a histogram of the time separations between the consecutive 
bursts. The result is shown in Fig. [2] The expected average separation is 64 days. The observed 
distribution roughly corresponds to the Poisson distribution of exponential decrease, as expected 
if the bursts are uncorrected in time. However, traces of enhanced small separations as well as 
big once may hint toward departure from truly random arrival. 

We also attempted a study of the arrival times separately for the bursts located in the Galactic 
Anticenter region (VSB1) and bursts located outside this region (VSB2). The number of bursts 
inside the region is enhanced at the beginning and seems to systematically decrease closer to the 
end of observation. The time distribution of bursts outside the Anticenter region seems to be more 
uniform. The^ 2 test of hypothesis of both distributions being constant in time give the fit quality 
of 11.0 and 13.0 per 8 d.o.f., for VSBi and VSB2, correspondingly. The Kolmogorov-Smirnov 
test gives the probability P = 0.567 of the VSBi distribution and whole VSB sample being the 
same. The distribution of the time separations for the whole VSB set and the bursts from the 
Galactic Anticenter region was compared using the Kolmogorov-Smirnov test. We multiplied 
the time separations for bursts from the Anticenter region by the ratio of the total number of 
bursts in the region to the total number of VSB. After such rescaling, the two distributions are 
consistent at the probability level of 0.976. 



2 http:/ ( 'http://f64.nsstc.nasa.gov/batse/grb/ 
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Figure 1: Distribution of the bursts arrival times for VSB. 
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Figure 2: Distribution of the bursts separation times for VSB. Expected distribution for random arrival times is marked 
with the dotted histogram. 
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Therefore we conclude that the bursts which come from the restricted region in the Galaxy, 
roughly towards the Galactic Anticenter, and the other VSB are not statistically different, and 
their arrival is statistically consistent with being uniform in time. 

2.2. Burst profiles 
2.2.1. Individual profiles 

We analyze the time profiles of the VSB from the event files for SB in the BATSE archival 
The requested data for 44 bursts of all VSB are available, including 19 bursts from the Anticenter 
region. We construct the lightcurve for each of the bursts using either all four channels, or 
selected channels. The individual profiles of bursts from the Galactic Anticenter region are shown 
in Fig. |3](all four channels), Fig. |4](channel 1) and Fig. [5] (channel 4), and bursts from all other 
regions of the sky are shown in Fig.[6](all four channels), Fig. [7] (channel 1) and Fig.[8](channel 
4). The bursts are labeled using their BATSE trigger number. 

All bursts show significant detections in the fourth channel. We determine the background in 
each channel separately and measure the peak countrate to the background ratio in 0.002 s bins. 
The ratio varies from 2.2 (for 6540) to 16.9 (for 5502). All VSB are hard, as discussed before by 
Cline et al. (2005) at the basis of KONUS data. 

Some of the bursts are clearly double-peaked (e.g. 2464, 3910, 6486, 6788), for one of the 
bursts (3412) the tte data starts after the peak, and for some bursts the data quality is so low that 
they cannot be used for the profile fitting analysis. 

Therefore, we perform the profile fitting only for 17 bursts, using all four energy channels. 
Additionally, we include the burst with the trigger number 0512 since it is very short although 
the formally given Tyo in BATSE catalog is somewhat higher than 0.1 s. 

We subtract the background for each of the bursts separately, and the background level is 
determined from the data points separated by more than 0. 15 s from the burst maximum, and we 
renormalize the bursts to have maximum at 1.0 for all events. 

We fit the normalized profiles using double exponential fit: 

/(f) = exp((f - f ) /fee) for t < t 

= exp(-(f - f )/ffaii) fort>t (1) 

and the three parameters, tn Se , f&n, and the shift to are fitted to the data. The time axis for each 
burst is defined in such a way that t — corresponds to the peak maximum (or to the first peak, 
if two equal count-rate peaks exist). We also made fits with fo = 0, i.e. with the change from 
rise to fall fixed at the maximum of the burst. However, better^ 2 is in general given by fits with 
to as a free parameter. In case of one of the bursts (BATSE trigger number 3904), a freedom 
of to shifts the peak considerably and the fit is problematic due to the fact that the background 
before the burst is much higher in this case than after the burst. Both fits are presented: selected 
bursts, normalized and background-subtracted, are shown in Figs. [9] and [10] for bursts from the 
Galactic Anticenter region and for those from the other parts of the sky, respectively. In the 
further analysis we use the fit with arbitrary to for all bursts for consistency. 

Almost all bursts show certain amount of asymmetry, with faster rise time and slower de- 
crease. The median of the rise timescale is 0.0040 s, the median of the fall timescale is longer, 
0.0110 s, and the median value of the ratio of these two timescales is 0.4. Therefore, the 



3 http://cossc. gsfc.nasa.gov/batse/batseburst/tte/ascii. html 
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Figure 3: Profiles of the VSB from the Galactic Anticenter region in channels 1-4. Count rate given for 0.002 s bins. 
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Figure 4: Profiles of the VSB from the Galactic Anticenter region in channel 1. Count rate given for 0.002 s bins. 
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Figure 5: Profiles of the VSB from the Galactic Anticenter region in channel 4. Count rate given for 0.002 s bins. 
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Figure 6: Profiles of the VSB from outside the Galactic Anticenter region in channels 1-4. Count rate given for 0.002 s 
bins. 
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Figure 7: Profiles of the VSB from outside the Galactic Anticenter region in channel 1. Count rate given for 0.002 s bins. 
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Figure 8: Profiles of the VSB from outside the Galactic Anticenter region in channel 4. Count rate given for 0.002 s bins. 
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Figure 9: Normalized background-subtracted profiles in 1 - 4 energy channels of the selected VSB from the Galactic 
Anticenter region in 0.002 s bins. The short-dashed line shows the double exponential fit to the burst profile with fixed 
tg = 0, while the long-dashed line shows the fit with to being a fitting parameter. 
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Figure 10: Normalized background-subtracted profiles in 1 - 4 energy channels of the selected VSB from outside of the 
Galactic Anticenter region in 0.002 s bins. The short-dashed line shows the double exponential fit to the burst profile 
with fixed to = 0, while the long-dashed line shows the fit with to being a fitting parameter. 
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Table 1 : Properties of single VSB with good S/N ratio 



BATSE trigger 

0512 
0207 
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6645 
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Figure 11: Histogram of the exponent ratio r = ffa.11 /^rise with f risc ar| d 'fall defined on the rising and falling part of the 
lightcurve of selected VSB 
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timescales on average are much shorter than found by McBreen et al. (2001) for short bursts. 
This is because they included also short bursts with Tgo > 0.1 s. The timescale ratio indi- 
cates larger asymmetry in selected VSB but this effect might be due to some contribution of 
double-peak events even in the selected bursts. For example, 2132 seems to have traces of a sec- 
ondary maximum at t as 0.015 s, and 3665 may actually be a multiple case. In Fig.QT|we show 
the distribution of the ratios of the two exponents in the two-exponent fit to the selected bursts 
(Equation Q]). The histogram has two peaks, one around r - ?faii/?rise=2.5 ratio and one around 
the ratio of r ~ 10. However, the effect of bimodality is not significant, taking into account the 
small number of bursts. 

We checked the burst profiles in all four energy channels for possible time delays. In some 
bursts the peak in the first channel was earlier than in the remaining three channels (a negative 
lag in trigger no. 0207 by 2 ms; in 0289 by 4 ms; in 0491 by 8 ms; in 2326 by 2 ms; in 3910 by 
4 ms), and in a few cases the peak in the channel 4th was delayed (a positive lag in trigger no. 
0432 by 12 ms; in 2159 by 2 ms; in 2332 by 8 ms; in 3384 by 8 ms; in 3803 by 6 ms; in 5536 by 
4 ms). In other bursts, the peak in channel 4 was advanced (a negative lag in trigger no. 2463 by 
2 ms; in 3037 by 4 ms; in 3810 by 2 ms;in 5502 by 2 ms; in 5620 by 4 ms; in 6547 by 2 ms; in 
7227 by 2 ms). In two bursts (2757 and 3073) two peaks in channel 4 are seen. In all other bursts 
the peaks in all energy channels coincided (with the accuracy of the 2 ms time bin) within the 
statistical error of one sigma. However the poor countrate prevents us from the detailed fitting of 
the profiles in separate energy channels and the quantitative measurement of the time delays was 
not possible, and the lags mentioned above are likely insignificant. 

We also compared the fluxes of the bursts within the Anticenter region and outside it. The 
bursts seen towards the Anticenter are somewhat brighter (with median countrate of 77 + 5 cts 
in 0.002 bins without background subtraction) than the remaining bursts (median of 66 + 6 cts), 
although the back ground itself determ ined close to VSB is actually slightly lower in the direction 



of the Anticenter dCline et all 120051 25.0 cts vs. 27.6 cts in 0.002 bins). However, since errors 



are large the difference is not highly significant. 

2.2.2. Composite profiles 

In order to overcome the problem of low S/N ratio in most of the bursts we also construct the 
composite burst profile. 

The background-subtracted lightcurve is later normalized by its maximum value and shifted, 
so the maximum is at t — 0. All lightcurves in all four channels are added, and the result divided 
by the number of lightcurves included. The resulting composite profile is shown in Fig. [T2] 

The profile shows significant asymmetry, characteristic also for longer bursts. Its shape is 
well represented by the exponential rise with the timescale of 0.0079 s, and approximately rep- 
resented by the exponential decay time of 0.0171 s. 

These timescales are much shorter than the average ris e and decay timescales of single pulses 



in a sample of SB (~ 0.036 s and ~ 5.5 s correspondingly, McBr een et al. L 2001) but fall within 



the broad distribution of these parameters i n the sample. Therefore, VSB show basically the 



same trend as seen in all gamma ray bursts (Band et all Il997t iNorris et all 120001) although the 



timescales involved are clearly much shorter. 

Better fit to the decay phase, however, is provided by the function 



/(f) = (irW (2) 
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Figure 12: Composite burst profiles for all VSB (black line), for bursts from the Galactic Anticenter region (red line) and 
for bursts from outside that region (blue line). The analytical fit (dashed line) is given by Eq.fJ Better fit for the decay 
part is provided by Ryde & Svensson function (Eq.ff] see text) 



with index n — 2.62 and the decay timescale t = 0.0319 s. iRvde and Svensson" (2002) analyzed 
22 bright pulses in long bursts using this description and they found a bi-modal distribution. 
Although most of bursts are well represented by n — 1 .0, there was a secondary enhancement 
around n = 3.0. It is therefore interesting to note that VSB may show some similarity to the latter 
subclass of pulses. 

If we choose only the energy range 50 - 300 keV (i.e. the channels 2+3) for the analysis the 
results are: timescale of the exponential rise 0.0093 s, exponential decay 0.0195 s, with again a 
much better fit to the decay part given by Eq.|2j with n = 2.44 and t = 0.0319 s. 

We searched for a correlation between the burst asymmetry and the burst hardness defined 
as a ratio of the peak countrate in the channel 4 to the total peak countrate. However, we do not 
detect a significant trend, perhaps due to the low quality of the data. 

We also calculated the composite burst profile in all four channels separately. We used the 
fixed time grid determined for the total profile (e.g. for combined 1-4 channels). In this way 
we allowed for advancement or a delay of the burst maximum with respect to the position of 
the maximum in the average profile from all four channels. The result is shown in Fig. [13] 
Burst peaks still coincide within the energy grid accuracy (0.002 s). We cannot increase the time 
resolution since the number of available photons is not large enough. The graph shows, however, 
that the profile in 4th energy channel is the narrowest. 

We also compare the shape of the composite pulse profile to spe cific prediction of the m odel 
of evaporation of a primordial black hole (for basic reviews, see e.g. ICarret allfl975l2005h . 

The time profile of the pulse produced due to the black hole evaporation can be derived from 
the temperature of the black body radiation, which in the Hawking process is inversely dependent 



16 




0.05 
time [s] 



Figure 13: Burst composite profiles for all VSB, for four energy channels separately, normalized to 1 for each energy 
channel: Red:channel 1, Greenxhannel 2, Bluexhannel 3, Black:channel 4. The narrowest profile is seen in channel 4. 
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Figure 14: Composite burst profile for all VSB (black line) with the analytical fits given by Eq. [6] The dotted line 
shows the fit with normalization calculated from the luminosity integral, while the dashed line shows the profile with 
normalization fitted to the background level. 
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on the black hole mass: 

T = — — — = 10 16 — !— [MeV]. (3) 
SnGkM M[g] 1 v ' 

The luminosity is then given by: 

dE d{Mc 2 ) A , 

L = = — = -o-T 4 4ttR 2 , (4) 

dt dt ' 

where R = 2GM/c 2 is the Schwarzschild radius. From the above relation, one can see that 
the luminosity is inversely proportional to the square of the black hole mass, M 2 . The mass 
dependence on time is found by integrating dM/dt from the time t < to fo = 0, when the black 
hole is completely evaporated. The lifetime of a black hole till evaporation is 

tocM 3 , (5) 

and the luminosity from Eq. |4]can be expressed as a function of time: 

L(t) oc cr 2/3 . (6) 

The luminosity profile goes to infinity at t — 0. Nevertheless, the integrated luminosity is fi- 
nite, so we can derive the normalization of the lightcurve based on the integral before the burst 
maximum being equal to that after the maximum. 

The fit of the function given by Eq. [6] to the composite profile is shown in Figure [14] Two 
curves represent the fits with (i) the normalization derived form the integrated luminosity (i.e. the 
integral of the analytic curve and the background-subtracted curve coincide) and (ii) the normal- 
ization arbitrarily fitted to the background level. The first method overpredicts significantly the 
luminosity of the burst before the peak. However, since the evaporation formula above does not 
give the time profile after the peak, the back ground determination may be biased by enhanced 
emission before as well as after the peak. This argument justifies the second approach, and in 
this case the fit to the rising part of the profile is good. 

The approach used above neglects the change of the spectral shape in time. It can be justified 
in the following way. The temperature of the black body emission during the final stage (~ 0.1 
s before complete evaporation) is of the order of 15 TeV, and at that stage a fireball is already 
formed. In other words, the evaporating plasma (photons, electron-positron pairs and heavier 
particles) is no longer optically thin, because a simple estimate of the total thickness for electron 
scattering indicates a depth of the order of 1 000. The following creation of electron-positron pairs 
is efficient and the final temperature of the radiation reaching an observer saturates at the energy 
treshold for this process. If this approximation is indeed justified then there is no temperature 
variation with time at the final stages and only the luminosity varies. Therefore the specific 
energy-dependent sensitivity of an instrument does not influence the shape of the lightcurve. 

Also, if we assume that indeed the emitted particles do not interact, the peak emission of 
photons is at very high energies, well above the instrument response. The detected photons 
would then come from the low energy tail, again perhaps weakly dependent on time apart from 
normalization. 

However, if the thermalizati on is not ful l , and for examp l e pion processes are im portant as 



discussed in several papers (e.g. iMaki et al.L 1996: Kapusta , 1999; Carr et all 120101) . or if the 
relativistic boosting of the fireball varies with time, there may be a spectral hardening effect 
and in this case our approach is too simplistic. However, no specific models of time-dependent 
spectral evolution of a fireball are available yet. 
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3. Discussion 



Very Short Gamma- ray Bursts were p roposed to originate as a result of evaporation of the 



primordial black holes (ICline et al.L 119991) . In this paper we carefully studied the spectral and 



temporal properties of VSB from BATSE with the aim to test whether they are likely to have a 
different origin than the other short gamma ray bursts. 

The burst duration for VSB is much shorter than the duration of short bursts, or of mixed 



samples with VSB and short bursts included (e.g. McB reen et al. Il2001h . The shortest rise time 
is only 5.3 ms. However, otherwise they seem to be the scaled-down analogues of short bursts. 
They are typically asymmetric, with faster rise than the decay time. Moreover, a number of VSB 
still consist of two or more subpulses, so the complex structure typical for longer bursts is visible 
even in this class. 

Significant fraction of VSB (21 out of 51) come from l/8th of the sky located roughly in 
the direction of the Galactic Anticenter region, extending from +90 to + 1 80 deg in galactic 



longitude and from -30 to + 30 deg in th e galactic latitude (ICline et aU 120031) . The significance 



of this anisotropy has been estimated by (Cline et al., 2001) to be 1.4 x 10 , and by Piotrowski 



etalQ (unpublished) to be a few times 10 using the factorial analysis. Bursts coming from the 
Anticenter region are slightly brighter but the effect is at the level of 2 sigma, and the distribution 
of the arrival times of both classes is statistically consistent with uniform. 

The enhancement of the number of VSB BATSE bursts in the direction of the Galactic Anti- 
center is clearly a puzzle. On the other hand, there are several interesting observations indicating 
the exceptional properties of that region. Roughly at that direction there is one of the richest 
known regions of star formation in the Galaxy, with over 800 HII regions, ma ny Wolf-Rayet 



stars, including the richest Cygnus OB2 association (see iKraemer et all 2010, for Cygnus-X 



Spicer Legacy Survey). The COMPTEL Al 26 survey indicated a cle ar enhancemen t in the An 



ticenter region also likely related to nucleosynthesis activity (see e.g. lDiehUl2001l) . An excess 



of the TeV cosmic rays from t he general d i rection of the heliotai l, also close to the Gala ctic An 



ticenter, has been reported bv lAbdo et alJ (12008b . However, (ISalvati and Sacco L 12008) argued 



against the heliotail origin. Similar claims of such cosmic ray anisotropy has also been done in 
the past, but wi t h the argument in favor of the connection with heliotail based on annual trend 



dFujimoto et all 1200 lb . Interesting possibility of correlations between the VSB and the regions 



of enhanced correlation betwe en the microwave back ground and the diffuse gamma rays above 
1 GeV has been considered bv ljedrzeiczak et alJ (120071) . 

For the energy distribution, VSB are much harder than short bursts (ICline et alll2005l) . How- 



ever, this effect can be considered as a continuation of a more general trend of bursts being harder 



when shorter, including the long-soft GRBs (see e.g. IZhang et all 12009b . We did not find any 



significant time delays between the peak emission in the four energy channels and we can only 
put an upper limit of 2 ms for such a delay. It is interesting to note that studies of long bursts 
showed that there i s an anticorrelation between the lags and the total burst isotropic luminosity 
20001) : 

\-1.15 



(Norrisetal. 



1.3( 



0.01* 



[10 53 ergs' 1 ] 



(7) 



It was interpreted by Salm onson and Gala ma (2002) as due to the variation in jet opening angles. 
They also suggested that measurements of variability or of spectral lags can then be used as a 
crude distance indicator for the GRBs. However, in the case of short bursts, both positive and 



4 http://www.lip.pt/events/2006/ecrs/proc/ecrs06-s5-10.pdf 
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negative lags were reported dGupta et al. , 2002), usually on the order of 50 ms. Therefore, the 
absence of measurable lags in combined profiles of VSB may not be surprising. 

The distribution of the GRBs in space shows that the cosmological effects are playing the im- 
portant role. The mean V/ V max value for the sample of long bursts was 0.282, while for the shor t 
GRBs it was somewhat larger, equal t o 0.390 ([Guetta and Piran , 2005 ; GuettaL and Piran C2006 ). 
In this respect, the value estimated by Cline et al.l (1 1999b to be 0.52. and indicating the euclidean 
distribution of the VSBs, seems to be a continuation of this trend. The shorter the bursts, the 
more locally they are observed. Based on our present study, is still not clear whether this might 
be the effect of a yet another different progenitor type, or rather a selection effect. 

3.1. VSB from SWIFT 

The number of very short events (T90 < 0. 1 s) detected by the SWIFT satellite is small so they 
are not suitable for statistical analysis. Nevertheless, these high quality data may shed some light 
onto the nature of such events. In Table [2] we list all the events with < 0.1 s, and in addition 
two more with Tgo just above 0.1. Out of those 10 sources, none is located in the Anticenter 
region found overpopulated wi th VSB in the BAT SE data. Three bursts now populate another of 
the 8 sky regions discussed bv lCline et al.1 (120011) . located between and 90 deg of the Galactic 
longitude, and -30 and +30 deg of the Galactic latitude. The probability of such coincidence is 
0.09 which cannot be considered as highly significant departure from uniform distribution. 

Candidate afterglow/host galaxies were claimed for some of those bursts but fading of the op- 
tical counterp art was never detect ed for any of those sources. For 100206A an optical counterpart 
was found by Levan et al. ( 2010l) . The spectroscopy obtained with the Keck telescope rev ealed 
the presence of a galaxy at the indicated location, with the redshift 0.41 ( Cenko et al.Ll2010h . and 
the projected offset of the burst from the galaxy center was ~ 35 kpc. The burst GRB 060502B 
has been found to coincide with z = 0.287 galaxy, and the implied offset from the center was 



73 + 19 kpc (IB loom et al.L 120071) . The event 050509B was very bright in X-rays and well local- 



ized, but again no fading optical/near infrared/millimeter counterpart was found. The source has 
been tentatively connected with z = 0.225 galaxy, a member of a galaxy cluster (IGehrels et al 



2005; iHiorth et al.Ll2005tlCastro-Tirado et alll2005l) . The projection distance in this case is ~ 33 
kpc. However, the X-ray afterglow measured by Swift XRT in this source start i ng from 62 s 



after the gamma burst was detected up to 300 s after the burst (IRol et al.L |2005; Bloo m et al 



20051) . and the X-ray luminosity was decreasing exponentially, with an index of ~ 1.3, implying 
the decay time below 60 s. The relatively long duration of this signal in comparison with the 
prompt gamma-ray emissi on may indeed support some (relativistic?) ejection of material which 
later undergoes cooling. iLee at al.l Q2005) found this timescale roughly consistent with the binary 
merger scenario. 

The comprehensive information on each of the bursts can be found through the internet 
for each of the bursts. The burst profiles viewed in 2 ms bins used for BATSE data analysis are 
usually single, although some substructure may be present in GRB 050925, GRB 0708 10B and 
GRB 0904 17 A. 



3.2. Possible mechanism of emission 

The tentative continuity of the burst properties from long bursts through short to VSB does 
not imply the same mechanism behind those events: long bursts are due to hypernova phe- 



5 http://www.swift. ac.uk 
6 http://gcn. gsfc.nasa.gov 
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Table 2: Properties of VSB from Swift 



name 


t90 

[s] 


gal. longitude 


galactic latitude 


redshift 


050509B 


0.073 


182.90828684 


86.15594183 


0.225 ? 


050925 


0.070 


72.30514279 


-0.07005621 




05 11 05 A 


0.093" 


59.86196122 


28.78790172 




060502B 


0.131 


81.59665275 


23.53577034 


0.287 ? 


070209 


0.090 


259.76518949 


-57.02797996 




07081 OB 


0.080 


116.37637070 


-53.85090173 




070923 


0.05 


295.85521498 


24.12470763 




0904 17A 


0.072 


173.46586705 


-61.03545070 




090515 


0.036 


232.46491111 


60.44320773 




100206A 


0.12 


166.88023385 


-37.72547634 


0.41 ? 


a t 

90 




0.028 


according 


to ( Cummings et al. 



2005), 



other 



burst parameters according 



to 



on-line SWIFT catalog 



http : //heasarc . gsf c . nasa.gov/docs/ swif t/archive/grb_table/ index . php. 



nomenon, short burst are due to mergers of compact objects in binary systems (neutron stars or 
black holes), and the origin of the VSB is still waiting for explanation. Two possibilities are 
natural to consider: a binary merger but with (for some reason) exceptionally short timescale and 
evaporation of primordial black holes. We consider those two scenarios below. 

3.2.1. binary merger 

The merger of binar y syste ms NS-NS and NS-BH were modeled by a number of authors 
(for the review, see e.g. Faberl l2009h . The expected duration of the BH-NS merger events is 



frequently somewhat longer than the shortest timescale 7*90 of a VSB (in our sample, the shortest 
VSB is 5.3 ms, see Table [Q in the SWIFT data the shortest burst found, GRB 090515, had T 90 
of 36 ms). 

The Newtonian simulations of the BH-NS merger (ILee and Kluzniak , 1999) were run for the 



mass ratio q - Mns/A^bh = 0.1 - 1.0, and showed that after about 30 ms the neutron star was 
disrupted and formed a torus of a mass of few tenths of the solar mass. The accretion rate at that 
moment was on the order of 2-6 solar masses per second, and therefore the expected lifetime of 
the torus feeding the black ho le in during the GRB event was T accr = M^/M = 40-60 ms. 
The relativistic simulations of Seti awan et al. ( 2006b . performed for a range of black hole spin 



parameters, followed the evolution of such torus with a mass of 0.01-0.2 M ,. They showed that 
in the end of the run, after 40-70 ms the mass of the torus is about half of the initial mass and the 
accretion rate drops to about 0.2-0.6 M /s. Therefore the lifetime of the central engine will be 
much longer than that. 

In merger event the duration of the phenomenon in the case of neutron star - black hole 
merger is fixed by the distance between the two at the phase of the neutron star tidal disruption 
and later by the viscosity oper ating within the accretion disk. 

In numerical simulation of IShibata and Taniguchi I d2008l) using T-law equation of state with 
r = 2 the disruption of the neutron star happened close to the innermost stable circular orbit. The 
disruption process lasts for about 5 ms and after that time the unswallowed remnants of the star 
form a torus of the size of about 50 km. The further stage of the accretion onto a black hole lasts 
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about 100 ms for a 3 M Q black hole, according to Janiu k et al. but those computations 



adopted 50 Rs for the size of the outer edge of the disk. If those computations are rescaled to 6 
Rs the duration of the accretion phase may possibly be considerably shorter, consistent with the 
decay timescale of the shortest burst. 

Short duration of the burst would then be provided by a small value of the black hole mass 
in the binary - larger black hole than 3 M would give too long burst duration. Thus in a natural 
way VSB would be connected with small mass black holes, would be fainter in the absolute 
luminosity and thus statistically less distant than short bursts. This last effect might explain the 
lack of cosmological effects in their V/V max distributio n. 



On the other hand, computations by iMetzger et al.l (2008) suggest that the duration of the 



NS-BH merger is longer if the disk mass (and in consequence, the total luminosity) is smaller, 
and for small disk masses (0.03 M Q ) the event lasts above 0.1 s, up to 1 s. If so, the VSB, likely 
to be intrinsically faint, are not consistent with VSB properties. 

Another merger scenario is possible if two neutron stars merge, and the total mass of the 
system is small. In this case it is probable that a differentially rotating hypermassive neutron 
star forms and initially is a stable configuration against the gravitational collapse. Later a de- 
laye d gravitational colla pse is in duced by the energy dissipation through gravitational waves, 



e.g. Shibata et al. ( 20061) : see also lOechslin et al" ( 2007 ). The collapse again leads to formation 



of hot, magnetized torus surrounding a rotating BH with the outer radius of high density region 
much smaller than 20 km, and in this case the expected lifetime of the torus in only ~ 1 ms. The 



detail ed predictions likely depend on the assumptions about the magnetic field (see e.g. lLiu et al 



2008) but in general the agreement of the predicted timescales and observed properties of VSB 



are in this case satisfactory. 

3.2.2. evaporation of primordial black holes 

The primordial black holes discussed here as a possible origin of the very short bursts should 
have masses in the range of 10 9 - 10 14 g, depending on time of their evaporation. The constraints 
for the abundance of the holes contribution to the dark matter density in the Universe, depending 



on their masses was recently discussed bv lLacki and Beacom I (12009) in a broad range of masses. 



The observed gamma-ray background limits the contrib ution of PBH to the critica l density of 
the Universe to Qpbh < 10~ 8 , whi ch was already s et by iPage and Hawking and simi- 



lar results are obtained nowadays (ICarr et all 12010). The limits for the somewhat larger PBH 
masses, obtained from the Milky Way neutrino, electron-positron and gamma ray backgrounds 
are constraining the total abundance of PBHs in the Universe and their contrib ution to the dark 



matter density, although the limits are less stringent dLacki and Beacom 1 120091 Qpbh < 10 for 
PBH masses between 2 x 10 15 g to 1000 M Q ). Abramowicz et al. (2009) studied the possible 
observational signatures of PBHs in the mass range of 10 16 - 10 26 g via their collisions with stars 
in our Galaxy, and found that the X-ray signals produced in such events would be too weak or 
too rare to be observable. 

The indentification of the VSB with the evaporating black holes does not violate those con- 
straints. If the signal is indeed due to evaporating black holes, the typical distance to the source 
can be estimated in the following way. The initial mass of the black hole evaporating now is of 
the order of 5 x 10 14 g, depending on the details of the evaporation efficiency at the later (hadron) 
stages. Nevertheless, most of the mass disappears during the initial slow evaporation. The mass 
of the black hole at the final stages, 0.1 s before the hole disappearance, is only about 6 x 10 8 g. 
If its energy is finally fully converted to photons and they reach BATSE detectors as ~ 100 keV 
the registration of a thousand photons would place the source at a distance of ~ 2 pc. Typical 
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distance estimates were a factor of 100 larger (see e.g. lChne et all 120051) since they were based 
on the initial PBH mass of 5 x 10 14 g. The numbers of VSB in the BATSE catalog, roughly 10 
per year, thus imply the average rate of 0.3 gamma ray flashes per cubic pc. We assume an initial 
steep power law distribution of the black hole masses (with a power law index larger than 2), 
and we take into account the mass evolution in time due to evaporation. The number of observed 
flashes is obtained from integration over the mass distribution from the minimum mass evaporat- 
ing now (6 x 10 8 g at the final 0.1 s) to the mass which will evaporate after a year (4 X 10 u g). 
Comparing this theoretical rate with the observed number of bursts we derive the normalization 
constant of the mass distribution. Integrating over the whole mass distribution, we obtain mass 
density of PBH equal 7 x 10~ 39 g cm 3 , i.e. Q. PBH ~ 10~ 9 . 

One of the expected properties of the black hole evaporation is the uniformity of the phe- 
nomenon. Black hole looses its angular momentum during the initial slow stage of the evap- 
oration (before the loss of 10 per cent of the mass). Thus the only parameter describing the 
evaporating black hole is its mass. If the number of particles released in a considered time unit is 
large, every event should look identically. The details of the emission process and the expected 
final radiation spectrum from t he fireba l l-type last stages have been considered in a number of 
papers (for the review see e.g. Bugaev , 2009). Although the predictions differ with respect to 
the energy peak of the spectrum as well as the spectral slope in the gamma-ray range, none of the 
papers predict any departure from isotropy of the emission, and consequently, from the identical 
look of all events. 

Our sample of VSB shows a variety of shapes, thus the events are not identical. This means 
that either the primordial black hole evaporation origin is ruled out, or the evaporation process is 
in need of a considerable revision. 

Such a revision is likely to come. As argued bv lMacGibbon et a l. (2008), the interaction be- 
tween the particles emitted by e vaporating black hole is probably not as strong, as co nsidered in a 



numb er of previous papers (e.g. lKapus ta , 2001; Daghigh and Kapusta , 2002); see Bu gaev et al 
(2008) for a review. If the particles finally escaping the gravity field of evaporating black hole 
move radially, since only almost purely radial velocity allows to leave the black hole vicinity, and 
the evaporation process proceeds through quark jets, there is indeed no thermalization involved 
at the beginning of the process and local fluctuations (i.e. in emission in various directions) may 
be large. The total energy of a VSB registered by an instrument (BATSE, in our case) is less 
than 1 erg per event, and thus all the registered emission of a given burst is likely to come from 
a single jet-like event close to a black hole horizon. If this is the case, then the events do not 
have to be exactly identical, prone to the initial quark ejection as well as tiny differences with 
respect to the inclina tion of a specific jet. H owever, they nevertheless should be short and hard, 
as already argued bv lCline and Hong I d 1992b . thus consistent with the VSB properties. 

We fitted to the composite VSB data the time profile of the optically thin fireball emitted 
in the final stage of the black hole evaporation. The rise timescale is in this case given by the 
analytic function, L(t) oc r 2 ^, while the normalization is fitted to the gamma ray background 
level. Such a model well represented the shape of the composite profile. After the maximum, 
the fireball will have different time characteristics, since we can have the hadronization with jet 
production, direct gamma and graviton production, as well as a direct lepton production, and 
large fluctuations in this process may explain the shape variety. This the black hole evaporation 
model is also consistent with the data for VSB. 

The tentative determination of the redshift for three VSB detected by SWIFT does not auto- 
matically rule out the evaporation scenario, because the connection with the candidate galaxies is 
only claimed by the proximity of the objects on the sky and may be a projection effect. However, 
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a long X-ray fading tail of the GRB 050509B discussed in Sect. 13.11 may indeed be difficult to 
explain within the evaporation scenario. 

4. Conclusions 

We studied the sample of the BATSE gamma ray bursts, limited to the special group of VSB 
(Tgo < 0.1 s). These bursts seem to originate in a special region in our Galaxy, however we 
found that their arrival is statistically consistent with being uniform in time. The time profiles of 
these bursts are asymmetric, with one or two peaks. We found no correlation between the burst 
asymmetry and its hardness, as well as no clear evidence for the time lags between the hard and 
soft photons. 

We found that the burst properties are most likely to be consistent with two origin mecha- 
nisms. The first plausible mechanism is the binary NS-NS merger with low total masses passing 
through a phase of hypermassive neutron star. The second one is the evaporation of a primordial 
black hole in the scenario of no photosphere formation. 
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